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Angular  intensity  correlations  in  the  double 
passage  of  waves  through  a  random  phase  screen 


Hector  M.  Escamilla,  Eugenio  R.  Mendez,  and  David  F.  Hotz 


The  problem  of  light  scattering  in  folded-path  or  double-passage  configurations  is  studied  theoretically. 
Assuming  as  the  random  medium  a  deep  phase  screen  that  introduces  Gauscian  distributed  phase 
fluctuations,  we  study  the  motion  of  the  speckle  as  the  source  is  moved.  Some  attention  is  also  given  to 
the  phenomenon  of  backscattering  enhancement.  Our  analysis  is  based  on  a  novel  expression  for  the 
complex  amplitude  that  has  a  simple  physical  interpretation.  For  simplicity,  only  the  one-dimensional 
case  is  considered,  but  an  extension  of  the  analysis  to  two-dimensional  screens  is  not  difficult. 

Using  the  factorization  properties  of  the  moments  of  a  complex  Gaussian  process,  we  are  able  to  derive 
analytical  expressions  for  the  mean  intensity  and  the  intensity  correlation  of  the  backscattered 
radiation.  We  find  that,  in  most  cases,  the  speckle  field  decorrelates  rapidly  as  one  moves  the  angle  of 
incidence  and  shifts  toward  the  direction  of  specular  reflection  with  a  rate  of  motion  that  is  different  from 
that  of  the  angle  of  incidence.  We  also  find  conditions  under  which,  when  the  angle  of  incidence  is 
modified,  the  speckle  pattern  produced  in  the  region  of  observation  tracks  the  backscattering  direction. 

Key  words:  Enhanced  backscattering,  do  le-passage  analysis,  speckle  motion,  multiple  scattering, 
rimdom  phase  screen. 


1.  Introduction 

The  problem  of  light  scattering  in  folded-path  or 
double-passage  configurations  has  been  studied  in  the 
last  three  decades  by  several  authors.  The  system 
is  known  to  exhibit  the  phenomenon  of  backscatter¬ 
ing  enhancement.  This  phenomenon  can  be  ex¬ 
plained  by  considering  the  coherent  addition  of  waves 
that  follow  identical  paths  in  opposite  directions  on 
interaction  with  the  random  medium.  Most  of  these 
studies  have  dealt  with  the  propagation  through  an 
extended  random  medium,  and  Ref.  6  provides  sev¬ 
eral  references  for  this  case.  More  recently,  the 
scattering  system  consisting  of  a  plane  mirror  in  the 
diffraction  field  of  a  deep  random  phase  screen  has 
also  been  studied  by  several  authors.'*  ®-’"'*’  To  our 
knowledge,  however,  only  the  mean  intensity  and  the 
first-order  spatial  coherence  function  of  the  scattered 
light  have  been  studied. 

In  this  paper  we  study  theoretically  the  deep  phase 
screen  double-passage  configuration  [see  Fig.  1(a)] 
with  particular  interest  in  the  correlation  of  the 
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intensity  fluctuations  that  describes  the  motion  of  the 
speckle  as  the  source  is  moved,  but  some  attention  is 
also  given  to  the  angular  distribution  of  the  mean 
intensity  around  the  backscattering  direction.  Our 
anedysis  is  based  on  am  expression  for  the  scattered 
complex  amplitude  that  can  be  obtained  by  manipulat¬ 
ing  the  diffraction  integrals  that  describe  the  propaga¬ 
tion  to  the  observation  point.  Assuming  that  a 
Gaussian  speckle  pattern”  is  formed  in  the  plane  of 
the  mirror  and  that  the  speckle  grain  on  the  mirror  is 
much  smaller  than  the  mirror  aperture,  we  derive 
exact  (within  the  model)  analytical  expressions  for 
the  mean  intensity  and  the  intensity  correlation  of 
the  backscattered  radiation  by  using  the  factorization 
properties  of  the  moments  of  a  complex  Gaussian 
process.  In  this  paper  we  deal  with  a  one-dimen¬ 
sional  random  phase  screen  that  introduces  smoothly 
varying  Gaussian-distributed  phase  fluctuations,  but 
the  calculations  can  be  easily  extended  to  the  case  of 
two-dimensional  phase  screens. 

We  find  that,  under  certain  circumstances,  when 
the  angle  of  incidence  is  modified,  the  speckle  pattern 
produced  in  the  region  of  observation  tracks  the 
backscattering  direction.  This  is  in  contrast  with 
the  well-  known  laws  of  speckle  motion  with  random 
surfaces*’-'®  and  with  the  type  of  motion  found  for 
volume  scatterers.'^'’  Another  interesting  feature 
present  in  our  calculations  is  that,  in  some  cases,  the 
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Fig.  1.  Scattering  geometry:  (a)  backscattering  configuration, 
(b)  equivalent  unfoHed  geometry  with  two  identical  screens. 


speckle  pattern  is  nearly  symmetric  about  the  back- 
scattering  direction. 

Using  these  results,  we  propose  a  modification  of 
the  original  scattering  geometry.  With  this  new 
configuration,  the  backscattering  enhancement,  the 
symmetry  of  the  speckle  pattern,  and  the  tracking  of 
the  specUe  pattern  of  the  backscattering  direction 
are  present  under  less  restrictive  conditions. 


2.  Formulation  of  the  Problem 

The  double-passage  configuration  is  shown  in  Fig. 
1(a).  A  plane  mirror  lies  in  the  Fresnel  region  of  the 
phase  screen.  Finite  ap)ertures  eu-e  located  in  the 
plane  of  the  mirror  and  of  the  screen.  The  angles  of 
incidence  0j  and  of  scattering  9,  su'e  taken  as  positive 
in  the  sense  indicated  by  the  arrows  in  the  figure. 
It  is  clear  that  the  unfolded  geometry  of  Fig.  1(b)  is 
equivalent  to  that  shown  in  Fig.  1(a).  The  aperture 
Wi  refers  to  the  width  of  the  incoming  beam. 


A.  Scattered  Complex  Amplitude 

With  reference  to  Fig.  1(b),  we  see  that  a  nearly 
collimated  monochromatic  Gaussian  beam  of  1/e 
amplitude  radius  Wi  is  incident  upon  the  first  phase 
screen,  which  is  located  on  the  plane  Xi,  at  an  angle  6, 
with  respect  to  the  z  axis.  This  incident  complex 
amplitude  Ai(jCi)  may  we  written  as 


A/xi)  =  Ko  exp(-ikxi  sin  0,)exp|-  -  (1) 

where  Kq  is  a  complex  constant  and  k  =  2tt/K,\  being 
the  wavelength  of  the  incident  light.  The  complex 
amplitude  of  the  scattered  light  in  the  plane  of  the 
mirror,  which  is  denoted  by  A(x,  0i),  can  be  found  by 
using  the  Kirchhoff  diffraction  theory  in  the  Fresnel 


approximation.  We  have 

A(x,  0,) 

r 

ik 

=  K,  I  A,(xi)exp|i4)(xi)]exp 

^(x,-x)2 

where  (J)(xi)  is  the  random  phase  introduced  by  the 
screen,  d  is  the  distance  from  the  screen  to  the 
mirror,  and  Ki  is  another  complex  constant.  The 
aperture  function  in  the  mirror  plane  is  given  by 

rn(x)  =  exp|- -  (3) 

where  W„.  is  the  aperture  radius  at  the  1  e  f)oint. 
The  complex  amplitude  A(x2),  which  is  incident  upon 
the  second  phase  screen  located  in  the  plane  X2,  can  be 
written  as 


A{X2) 


ik 

J  A(x,  0,)ro(x)exp 

^(X2-x)2 

dx,  (4) 


where  K2  is  a  complex  constant,  and  we  have  again 
made  use  of  the  Fresnel  approximation.  At  an  obser¬ 
vation  point  in  the  far  field  of  the  second  screen,  the 
scattered  complex  amplitude  can  be  written  as 

A(0„  ej  =  kJ  A(x2)exp|- 

X  exp[i<t>(x2)jexpl  -ikx2  sin  djdx2,  (5) 

where  0s  denotes  the  direction  of  observation  and  K  is 
another  complex  constant;  we  have  included  a  Gaus¬ 
sian  aperture  function  of  radius  W2  in  the  plane  of  the 
second  screen.  From  Eqs.  (3)-(5)  we  may  write, 
after  interchanging  the  orders  of  integration. 


A(0, 
where 
A(x 


',,63)  =  a:  I 


expi  -  ^ 


A(x,  0,)A(x,  0s)dx,  (6) 


,ej  =  K2j 


ik 


X  exp  ^(x2-x)^|exp(-jA!X2sin0s)dx2.  (7) 


The  quantity  A(x,  0s)  is  the  complex  amplitude  that 
a  point  source  located  at  the  observation  point  in  the 
far  field  would  produce  in  the  plane  of  the  mirror. 
Substitution  of  Eq.  (1)  into  Eq.  (2)  shows  that,  except 
for  an  irrelevant  constant  factor,  Eqs.  (2)  and  (7)  have 
exactly  the  same  form.  Our  approach  is  based  on 
Eq.  (6)  and  a  few  additional  assumptions  that  are 
described  in  Subsection  2.B.  Equation  (6)  has  a 
simple  physical  interpretation:  mathematically  it  is 
the  product  of  two  speckle  patterns  integrated  over 


20  May  1993  /  Vol.  32,  No.  15  /  APPLIED  OPTICS 


2735 


the  mirror  aperture.  This  expression  has  some  close 
parallels  with  the  equation  describing  image  forma¬ 
tion  in  confocal  scanning  microscopy. 

B.  Correlation  of  the  Scattered  Intensity 
At  this  point  it  is  necessary  to  introduce  a  statistical 
model  for  the  phase  fluctuations  <t)(x).  We  assume 
that  ({)(x)  is  a  stationary  zero-meem  Gaussian  random 
process  with  a  Gaussian  correlation  coefficient 
p(jc  -  jc'); 


we  can  write 

e,;  e,',  e;)  =  c,(e„  e,;  e/,  e;)  +  Caie,,  e,;  e,',  e;), 

(12) 

where 

Ci(8„e,;e.',0;) 


p(x  -  x') 


(<t>(x)(l)(x')) 


/ 


>  (8) 


and 


X  (A(x,  es)A*(x',0s'))dxdx’,  (13) 


where  is  the  variance  of  the  random  phase  fluctua¬ 
tions  and  ^  is  the  characteristic  correlaton  length  of 
the  process. 

Now,  for  a  deep  phase  screen,  <t^  »  1,  and  if,  in 
addition,  and  (kt^/2d)  «s:  1,  then  the 

complex  amplitude  A(x,  0,)  corresponds  to  that  of  a 
fully  developed  speckle  pattern."*’^'  That  is,  A(x,  8,) 
is  a  zero-mean  circular  complex  Gaussian  random 
variable  (CCGRV).“  The  condition  (ki^/2d)  «:  1 
ensures  that  the  mirror  is  far  from  the  focusing 
region  of  the  screen  where  strong  departures  from 
Gaussian  statistics  take  place.  Likewise,  if  W2^  » 
^2,  A(x,  0s)  is  also  a  zero-mean  CCGRV.  In  addition, 
if  W„  is  much  greater  than  the  speckle  size  in  the 
plane  of  the  mirror,  the  scattered  complex  amplitude 
A(0,  ,  0s)  is  also  a  zero-mean  CCGRV. 

Consider  now  the  intensity  /(Oj,  0,)  at  detection 
angle  0,,  which  is  due  to  an  incoming  beam  incident  at 
angle  0,  ,  and  the  intensity  7(0,  ',  0s')  at  detection  angle 
0,',  which  is  due  to  an  incoming  beam  at  incident 
angle  d/.  With  the  statements  made  above  regard¬ 
ing  the  statistics  of  A(0„  0,),  the  intensity  correlation 
C/(0„  0s;  0,',  0s')  =  (/(O.,  0s)^(0.'.  0s'))  can  be  written 
as^”; 

C,(0„  0s;  0,',  0s') 

=  (7(0.,  0s)){7(0,',  0s'))  +  |Cx(0„  0s;  0,',  0s') P,  (9) 


C2(0„8s;0,',0s') 

=  |7:pJ  j*  exp|-^^^;j^j(A(x,0,)A*(x',0s')) 

X  {A(x,ejA*(x',e/)}dxdx',  (14) 

The  last  two  expressions  involve  only  the  calcula¬ 
tion  of  second-order  moments  or  correlations  of  the 
held  amplitude,  and  are  much  simpler  to  calculate 
than  the  fourth-order  moment  of  Eq.  (11).  The 
amplitude  correlations  in  Eqs.  (13)  and  (14)  were 
evaluated  in  the  limit  »  1  by  using  a  steepest- 
descents  method  (see  Appendix  A).  Using  these 
results,  we  can  obtain  analytical  expressions  for 
Ci(0i,  0,;  0,',  0s')  and  C2(0„  0*;  0,',  0s')  without  further 
approximations.  The  result  involves  some  rather 
long  algebraic  expressions  that  are  difficult  to  reduce 
and  thus  are  given  in  Appendix  A.  A  simple  formula 
for  a  particular  case  is  given  in  Section  3. 

The  mean  intensity  is  then  calculated  from 

(7(0,,  0s))  =  C^(8i,  0,;  0„  8.),  (15) 

while  the  normalized  intensity  correlation  ■yi(8„  0,; 
0/,  0,')  is  given  by 

7.(0,,  0s;  0,',  0s')  =  1  +  1 7x(0„  0s;  0.',  0s')  P 


where 

C’a(0„  0s;  0,',  0s')  =  (A(0„  0s)A*(0.',  0s')>.  (10) 

Using  Eq.  (6),  we  see  that  the  problem  reduces  to  the 
calculation  of  the  function: 

Ca(0.,0s;0.',0s') 


X  (A(x,0,)A(x,0,)A*(x',0/)A*(x',0s'))dxdx'.  (11) 

Recalling  again  the  statistical  considerations  made 
above,  we  also  find  that  the  fourth-order  moment  in 
the  integrand  of  Eq.  ( 1 1)  can  be  expressed  as  a  sum  of 
products  of  first-order  correlation  functions.^®  Thus 


|C,4(e.,es;0.',0s')P 

C^(0,,0,;0i,0s)CA(0i',0s';0.',0s')' 

(16) 

It  is  pertinent  to  mention  that,  recently,  a  compact 
formula  for  the  mean  intensity  produced  by  the 
scattering  of  a  Gaussian  beam  by  the  double  passage 
through  a  deep  phase  screen  has  been  derived  that 
also  uses  the  factorization  properties  of  the  moments 
of  a  Gaussian  random  process.*® 

3.  Resutts 

We  now  turn  our  attention  to  the  study  of  the 
behavior  of  the  angular  distribution  of  the  mean 
intensity  around  the  backscattering  direction.  Fig¬ 
ure  2  shows  the  evolution  of  the  enhancement  factor 
Ef  as  the  distance  d  from  the  phase  screen  to  the 
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Fig.  2.  Enhancement  factor  Ef  as  a  function  of  the  normalized 
distance  d  L  between  the  mirror  and  the  phase  screen.  For  the 
numerical  calculations  VF2  and  were  assumed  to  be  infinite,  L  = 
6.206  mm  (see  text)  and  6,  =  5  mrad. 


mirror  increases  for  the  case  of  infinite  apertures  at 
the  mirror  and  the  second  screen.  We  define  the 
enhancement  factor  Ef  as  the  ratio  of  the  mean 
intensity  in  the  backscattering  direction  to  the  mean 
background  intensity.  In  Fig.  2,  the  circles  repre¬ 
sent  the  values  obtained  from  the  evaluation  of  the 
general  formulas  given  in  Appendix  A,  while  the  solid 
curve  is  a  plot  of  the  expression 


Ef  =  1  -f 


(17) 


‘  '  [1  +  (d/L)‘^Y  2 

where  the  normalization  distance  L  is  given  by 

L  =  (18) 


2a 


For  the  numerical  calculations  presented  we  have  set 
\  =  0.6328  nm,  a  =  8  rad,  4  =  5.0  p,m,  and  Wj  =  2 
mm.  So  the  paraimeter  L  was  fixed  at  the  value  L  - 
6.206  mm.  Equation  (17)  can  be  derived  from  the 
results  given  in  Ref.  4.  The  agreement  between  the 
two  plots  is  excellent.  It  is  seen  that  the  height  of 
the  peak  diminishes  as  the  screen  to  the  mirror 
distance  increases.  This  is  because,  as  this  distance 
increases,  more  of  the  scattered  light  reflected  by  the 
mirror  onto  the  screen  lies  outside  the  region  illumi¬ 
nated  in  the  first  passage.  It  is  only  the  region  of 
overlap  that  contributes  to  the  backscattering  en¬ 
hancement.  The  light  reflected  outside  this  region 
contributes  to  only  the  background  scattered  intensity. 

The  relative  enhancement  can  also  be  increased  by 
reducing  the  area  of  the  screen  illuminated  in  the 
second  passage  by  placing  a  physical  aperture  at 
either  the  mirror  or  at  the  screen.  This  is  illustrated 
in  Fig.  3  for  a  normalized  distance  d/L  =  6.5,  for 
which  expression  (17)  predicts  a  small  enhancement 
(see  Fig.  2).  In  Fig.  3(a)  the  aperture  radius  of  the 
mirror  has  been  set  to  one  half  of  the  radius  of  the 
incoming  Gaussian  beam.  A  high  value  is  obtained 
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Fig.  3.  Angular  distribution  of  mean  intensity  as  a  function  of  the 
scattering  angle  9,  for  the  cases  of  reduced  apertures  at  the  mirror 
or  phase  screen.  The  dots  indicate  the  direction  of  incidence  B,  =  5 
mrad.  Hered  =  6.5  L  and  (a)  W„  =  Wj  2  and  infinite,  (b)  W2  = 
VFi  and  W„  infinite. 
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for  the  relative  enhancement  because,  after  reflec¬ 
tion,  the  illuminated  region  of  the  second  screen 
almost  coincides  with  the  area  illuminated  in  the  first 
passage.  Figure  3(b),  on  the  other  hand,  shows  that 
the  maximum  theoretical  enhancement  factor  of  2 
(see  Refs.  21-23)  is  achieved  when  the  apertures  at 
the  two  screens  have  equal  sizes,  as  this  maximizes 
the  proportion  of  reversible  paths.  The  mirror,  in 
this  case,  has  been  assumed  to  be  infinite. 

In  Fig.  4  we  study  the  motion  of  the  speckle  pattern 
as  the  direction  of  the  incoming  beam  changes.  We 
proceed  by  first  fixing  the  original  angles  of  incidence 
0i  and  detection  0,.  We  then  choose  a  new  angle  of 
incidence  0,'  and  find  the  maximum  value  of  the 
correlation  7,(0,,  0,;  0,',  0,')  as  we  scan  the  second 
angle  of  observation  0/.  In  these  figures,  the  verti¬ 
cal  dashed-dotted  line  indicates  the  original  direction 
of  observation  0,  while  the  dotted  line  indicates  the 
original  direction  of  incidence  0,.  Figure  4(a)  is  a  plot 
of  the  maximum  value  of  the  intensity  correlation 
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Scattering  Angle  (mrad) 


Incidence  Angle  (mrad) 

Fig.  4.  Envelopes  of  the  intensity  correlation  functions.  The 
dashed-dotted  vertical  line  in  (a)  indicates  the  original  direction  of 
observation,  while  the  dotted  line  in  (b)  indicates  the  original 
direction  of  incidence.  Here  d  =  6.5  L,  =  W-i,  is  infinite, 
e,  =  0  mrad  and  0,  =  -1  mrad.  (a)  Maximum  values  of  the 
intensity  correlation  function  as  a  function  of  the  scattering  angle 
0,',  lb)  maximum  values  of  the  intensity  correlation  function  as  a 
function  of  the  angle  of  incidence  0,'  (circlesi  and  angular  distribu¬ 
tion  of  mean  intensity  with  0,'  =  0  mrad  as  a  ‘‘unction  of  the  angle 
of  incidence  6,’  (solid  curve). 


function  as  a  function  of  the  angle  of  observation  S,' 
at  which  the  maximum  occurs  for  each  new  angle  of 
incidence  0,  ’;  the  curve  represents,  then,  the  intensity 
correlation  envelope.  It  should  be  pointed  out,  how¬ 
ever,  that,  in  general,  the  ratio  of  the  angular  shift  of 
the  speckle  pattern  to  the  angular  shift  in  the  angle  of 
incidence  is  not  equal  to  1 .  This  ratio  depends  on  the 
geometry  and,  in  particular,  it  depends  on  the  dis¬ 
tance  d.  To  illustrate  the  fact  that  the  speckle  does 
not  move  at  the  same  rate  as  the  angle  of  incidence, 
we  show  in  Fig.  4(b)  a  plot  of  the  correlation  envelope 
shown  in  Fig.  4(a)  in  terms  of  the  direction  of 
incidence  0,'.  The  correlation  envelope  shown  in 
Fig.  4(b)  then  represents  (for  given  angles  0„  0,,  and 
0,')  the  maximum  correlation  obtained  by  scanning 


the  angle  of  scattering  0/  plotted  as  a  function  of  the 
angle  of  incidence  0,'.  It  can  be  seen  from  Fig.  4(a) 
that,  for  the  distance  d  =  6.5  L  employed  in  the 
figure,  the  full  width  of  the  correlation  envelope,  as  a 
function  of  the  scattering  angle  0/,  equals  approxi¬ 
mately  2  M.rad.  This  width  corresponds  to  an  inter¬ 
val  of  angles  of  incidence  0,'  of  approximately  3  mrad, 
as  can  be  seen  from  the  curve  plotted  in  Fig.  4(b). 
This  indicates  that,  in  the  situation  considered,  the 
speckle  pattern  is  practically  static,  and  decorrelates 
as  the  angle  of  incidence  is  changed.  The  sofid  curve 
in  Fig.  4(b)  represents  the  mean  intensity  for  0/  =  0 
mrad  (which  corresponds  to  the  original  angle  of 
incidence  0,)  plotted  as  a  function  of  the  angle  of 
incidence  0,’.  Invoking  reciprocity  arguments  we 
can  see  that  the  same  curve  would  be  obtained  by 
setting  0,  =  0  mrad  and  plotting  (/(0,,  0s))  as  a  function 
of  Oj.  The  solid  curve  in  Fig.  4(b)  represents,  then, 
the  backscattering  enhancement  peak,  emd  it  is  re¬ 
markable  that  the  two  curves  plotted  in  this  figure 
are  practically  the  same.  A  similar  result  holds  for 
scattering  from  dense  media,  whe’'“  it  has  been 
shown  that  the  shape  of  the  envelope  ol  the  <’orrela- 
tion  is  equal  to  the  square  of  the  backscattering 
enhancement  peak.^-*  We  also  mention  that  the  situ¬ 
ation  shown  in  Fig.  4(b)  does  not  always  hold.  For 
instance,  if  the  mirror  is  moved  closer  to  the  phase 
screen,  one  eventually  reaches  a  point  beyond  which 
the  shape  of  the  envelope  correlation  and  that  of  the 
backscattering  peak  are  different. 

Figure  5  shows  two  aspects  of  the  behavior  of  the 
intensity  correlation  function  as  the  direction  of 
incidence  changes  from  0,  to  0, ' .  The  geometry  of  the 
arrangement  is  the  same  as  in  Fig.  4,  but  we  have  now 
set  the  distance  d  =  0.8  L.  Also,  as  in  Fig.  4,  the 
vertical  dashed-dotted  line  in  the  upper  half  of  these 
figures  indicates  the  original  direction  of  observation 
05,  while  the  dotted  lines  indicate  the  original  direc¬ 
tion  of  incidence  0, .  The  long-dashed  vertical  lines  in 
the  bottom  hedves  of  the  figures  indicate  the  new 
direction  of  incidence  0,',  while  the  solid-vertical  lines 
indicate  the  new  direction  of  observation  O,'  for  which 
maximum  correlation  is  achieved.  In  Fig.  5(a),  the 
thick  solid  curve  corresponds  to  the  intensity  correla¬ 
tion  function,  while  the  thin  curve  shows  the  enve¬ 
lope  of  the  intensity  correlation  function.  For  the 
value  of  d/L  employed  in  Fig.  5(a),  the  envelope  is 
much  wider  than  that  in  Fig.  4(a).  We  also  see  from 
Fig.  5(a)  that,  as  the  direction  of  incidence  changes 
from  0,  =  0.0  mrad  to  0/  =  0.15  mrad,  the  peak  of  the 
intensity  correlation  function,  and  thus  the  speckle 
pattern,  moves  in  the  direction  of  the  specular  reflec¬ 
tion  from  the  original  angle  of  observation  6,  =  - 1 
mrad  to  the  position  0s'  =  - 1.04  mrad.  In  Fig.  5(b), 
the  original  directions  of  incidence  eind  observation 
are,  respectively,  0^  =  0.0  mrad  and  0,  =  -1.0  mrad. 
It  can  be  seen  that,  when  the  new  direction  of 
incidence  0,'  coincides  with  the  original  direction  of 
observation  0,,  the  peak  of  the  intensity  correlation 
function  appears  at  the  original  direction  of  incidence 
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Fig.  5.  Shift  of  the  intensity  correlation  peak  with  a  change  in  the 
angle  of  incidence.  The  abscissas  correspond  to  the  scattering 
angle  9*'.  The  dashed-dotted  vertical  line  indicates  the  original 
direction  of  observation,  while  the  dots  indicate  the  original 
direction  of  incidence.  The  long-dashed  vertical  line  indicates  the 
new  direction  of  incidence  while  the  solid  vertical  line  indicates  the 
angle  of  observation  for  which  maximum  correlation  occurs. 
Here  d  =  0.8  L,  Wi  =  W2,  W„  is  infinite,  0;  =  0  mrad,  and  9,  =  - 1 
mrad.  (a)  Intensity  correlation  function  for  0/  =  0.15  mrad  (thick 
curve'  and  envelope  of  the  intensity  correlation  function  (thin 
curve).  The  maximum  correlation  occurs  at  0,  =  -1.04  mrad. 
(b)  Intensity  correlation  function  for  0^'  =  0,. 


0;.  So,  when  source  and  detector  are  interchanged, 
the  corresponding  intensities  are  perfectly  correlated. 
This  result  shows  that  the  formulas  derived  in  Appen¬ 
dix  A  obey  reciprocity.  This  kind  of  correlation  was 
recently  discussed  for  volume  scatterers  and  has  been 
termed  the  “time-reversed  memory  effect.”*'^ 

In  Fig.  6  we  consider  the  case  of  a  small  aperture  at 
the  mirror,  much  smaller  than  the  cross  section  of  the 
incoming  beam.  For  these  figures  we  assumed  an 
infinitely  large  aperture  at  the  second  screen.  Fig¬ 
ure  6(a)  shows  the  spatied  intensity  autocorrelation 
function.  It  can  be  seen  that  the  speckle  pattern  is 
nearly  symmetrical  with  respect  to  the  back^ttering 


Scattering  Angle  (mrad) 


Scattering  Angle  (mrad) 

Fig.  6.  Intensity  correlation  functions  as  functions  of  the  scatter¬ 
ing  angle  0,'  for  a  small  mirror  aperture.  The  vertical  lines  have 
the  same  meaning  as  in  Fig.  5.  Here  d  =  L,W„  =  0.025  W],  W2  is 
infinite,  0;  =  0  mrad,  and  9,  =  -1  mrad.  (a)  Intensity  correlation 
for  9j'  =  0,  (spatial  autocorrelation),  (b)  intensity  correlation 
function  for  9/  =  1.0  mrad. 

direction.  This  is  because  [with  reference  to  Fig. 
1(b)]  the  light  scattered  from  the  first  screen  that 
reaches  the  second  screen  goes  through  a  small 
aperture  on  the  axis  of  the  system  so  that  an  inverted 
image  of  the  area  illuminated  on  the  first  screen  is 
projected  onto  the  second.  Figure  6(b)  shows  that 
when  the  incident  direction  changes  from  0,  =  0  to 
0,'  =  1  mrad,  the  correlation  peaks  move  in  a  direc¬ 
tion  opposite  to  that  of  the  specular  direction  and  at 
the  same  rate  as  the  incident  direction.  That  is,  in 
this  case,  the  speckle  pattern  tracks  the  backscatter- 
ing  direction. 

These  striking  results  are  due  to  the  fact  that,  for 
the  complex  amplitudes  A(x,  0j)  and  A(x,  0,)  [see  Eq. 
(6)],  the  phase  curvatures  within  the  mirror  aperture 
are  negligible.  This  su^sts  a  modification  of  the 
scattering  geometery  that  will  produce  these  effects 
under  less  restrictive  conditions.  The  new  geometiy 
is  shown  in  Fig.  7.  The  simple  formula  referred  to  in 
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Fig.  7.  Double-passage  configuration  with  the  phase  screen  at  the 
center  of  curvature  of  a  spherical  mirror. 


Subsection  2.B  was  derived  for  the  case  =  W2  =  W'’ 
and  for  this  specialized  geometry,  including  a  spheri¬ 
cal  mirror,  with  an  aperture  of  radius  W^,  situated  in 
the  feir  field  of  the  screen.  The  radius  of  curvature  of 
the  mirror  is  d  and  is  equal  to  the  distance  from  the 
screen  to  the  mirror. 

It  is  found  that  up  to  a  constant  factor,  the 
amplitude  correlation  is  given  by  the  expression 


7a(0..  0s;  0.'.  0s') 

1  (e>  +  o/)2  1  (6,-  +  0,-)2 


=  exp 
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dieted  by  this  equation  is  given  by 
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which  predicts  a  broad  envelope  with  a  narrow  peak 
(width  (ps)  in  the  backscattering  direction  (0^  =  0,). 

Figure  8  corresponds  to  the  geometiy  depicted  in 
Fig.  7,  in  which  the  phase  screen  is  at  the  center  of 
curvature  of  the  spherical  mirror.  Figure  8(a)  shows 
the  behavior  of  the  mean  intensity.  It  is  seen  that 
the  relative  value  of  the  backscattering  enhancement 
reaches  a  value  of  2.  The  reason  for  this  is  that,  in 
the  present  case,  an  image  of  the  region  of  the  screen 
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These  last  three  quantities  represent,  respectively, 
the  angular  widths  of  the  speckle,  of  the  speckle 
envelope,  and  of  the  mirror.  To  derive  expression 
(19)  we  assumed  that  (p^,  (pm  »  Ps- 
Also,  it  is  straightforward  to  show  that,  with 
approximations  consistent  with  those  employed  in 
the  derivation  of  Eq.  (19),  the  mean  intensity  pre- 


Scattering  Angle  (mrad) 

Fig.  8.  Mean  intensity  and  intensity  correlation  for  the  geometry 
shown  in  Fig.  7.  The  vertical  lines  have  the  same  meaning  as  in 
Fig.  5.  The  relevant  partuneters  are  d  =  80.0  L,  W„  =  5  Wi,  W2  is 
infinite,  and  0,  =  5  mrad.  (a)  Angular  distribution  of  mean 
intensity  as  a  function  of  0,,  (b)  intensity  correlation  plotted  as  a 
function  of  the  scattering  angle  0,',  with  0,  =  2  mrad  and  0/  =  6 
mrad. 
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illuminated  by  the  incoming  beam  is  formed  on  the 
screen  by  the  mirror.  This  image  is  inverted  and  is 
of  unit  magnification.  The  intensity  correlation  func¬ 
tion  in  Fig.  8(b)  shows  the  same  features  obcerved  in 
Fig.  6.  The  autocorrelation  function  shows  two 
peaks,  indicating  the  symmetry  of  the  speckle  pattern 
,  :th  respect  to  the  backscattering  direction.  Also, 
as  the  angle  of  incidence  of  the  incoming  beam  is 
changed,  the  speckle  pattern  tracks  the  backscatter¬ 
ing  direction  and  not  the  specular  direction,  as  would 
be  the  case  with  a  single-scattering  rough  surface. 
These  features  can  be  explained  from  Eq.  (19):  The 
term(6r  -  B,  )  in  this  equation  is  the  one  responsible 
for  the  tracking  of  the  backscattering  direction,  as  it 
predicts  a  peak  when  sin  6,'  =  sin  0.'  +  (sin  0,  - 
sin  0,).  On  the  other  hand,  the  term  (0,^  -  0,^)  is 
responsible  for  the  symmetry  of  the  speckle  pattern 
about  the  backscattering  direction  (i.e.,  0/  =  0,'),  as  it 
predicts  a  peak  when  sin  0,'  =  sin  0,'  -  (sin  0,  - 
sin  0.). 

4.  Summary  and  Conclusions 

We  have  studied  the  scattering  of  light  resulting  from 
double  passage  through  a  deep  random  phase  screen 
for  an  incoming  Gaussian  beam.  Our  analysis  is  for 
a  one-dimensional  random  phase  screen  but  the 
extension  to  two  dimensions  is  streiightforward. 
Assuming  that  the  complex  amplitude  of  the  speckle 
formed  on  the  mirror  follows  Gaussian  statistics,  and 
that  the  speckle  grain  on  the  mirror  is  much  smaller 
than  the  mirror  aperture,  we  have  derived  formulas 
for  the  mean  intensity  of  the  backscattered  light  and 
for  the  correlation  of  the  intensities  for  different 
directions  of  incidence  and  observation.  Using  the 
factorization  properties  of  the  moments  of  a  complex 
Gaussian  process,  we  have  derived  exact  (within  the 
model)  analytical  expressions  for  these  two  quantities. 
The  formulas  so  obtained  are  long  and  difficult  to 
interpret  but  we  have  studied  a  variety  of  double¬ 
passage  configurations  from  the  numerical  evedua- 
tion  of  these  erp’^essions. 

We  have  studied  the  backscattering  enhancement 
factor  as  well  as  the  motion  of  the  speckle  pattern  as 
the  angle  of  incidence  is  changed.  We  have  found 
that,  in  general,  the  rate  of  change  in  the  position  of 
the  speckle  pattern  is  not  the  same  as  that  of  the 
angle  of  incidence.  In  normal  circumstances,  the 
speckle  moves  in  the  direction  of  the  specular  reflec¬ 
tion  at  a  rate  slower  than  that  of  the  incident  beam 
and  decorrelates  with  em  envelope  described  approxi¬ 
mately  by  the  shape  of  the  backscattering  enhance¬ 
ment  pe^.  Conditions  have  also  been  found  under 
which  the  speckle  pattern  is  symmetric  about  the 
backscattering  direction  and  tracks  the  backscatter¬ 
ing  direction  as  the  source  is  moved.  The  conditions 
for  the  observation  of  this  effect  are,  however,  fairly 
restrictive.  The  new  geometry  proposed,  with  a 
spherical  mirror,  exhibits  enhanc^  backscattering, 
the  symmetry  of  the  speckle  pattern,  and  the  tracking 
of  the  backscattering  direction  under  less  restrictive 
conditions. 


Appendix  A 

Apart  from  some  constant  factors,  the  amplitude 
correlations  appearing  in  Eqs.  ( 13 )  and  ( 14 )  are  of  the 
general  form 

(A(x,0)A*(x',0')) 

ik 

X  exp  -  xf  -  [Xj  -  x)^] 

X  exp{-iA!(x,  sin0  -x^sin0')]dx.dxj.  (Al) 

The  steepest-descent  method'“  “  that  we  employed 
for  the  calculation  of  these  expressions  basic^ly 
consists  of  making  the  approximation 
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in  the  evaluation  of  the  amplitude  correlations  given 
by  Eq.  (A.1).  Using  these  results,  we  find  that  Eq. 
(13)  gives 
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Similarly,  it  is  found  that  Eq.  1 14)  gives  the  follow¬ 
ing  result; 
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The  expressions  corresponding  to  the  geometiy 
depicted  in  Fig.  7  are  u’otained  by  dropping  the 
imaginary  terms  inside  parentheses  in  Eqs.  (All)  and 
)A21).  Equation  119)  was  obtained  by  dropping,  in 
addition,  the  imaginary  terms  in  Eqs.  (A12i  and  (A13i 
and  by  putting  IFi  =  VF2  =  VF. 

This  work  has  been  supported  by  an  internal  grant 
from  the  Centro  de  Investigacion  Cientifica  y  de 
Educacion  Superior  de  Ensenada  and  by  the  U  S. 
Army  Research  Office  under  grant  DAAL03-89-C- 
0036. 

References 

1.  N  G  Denisov,  “On  wave  scattering  under  full  reflection 
conditions,"  R,uiioph>s.  Quantum  Electron  7,  378-380  1964 : 

2.  A.  G  Vinogradov.  Yu.  A.  Kravstov.  and  V  I  Tatarskii.  The 
effect  of  intensification  of  backscattenng  by  bodies  placed  m  a 
medium  with  random  inhomogeneities. "  Radiophys  Quantum 
Electron.  16,  818-823  1 1973 1. 

3.  Yu.  A.  Kravstov  and  A.  1,  Saichev,  "Effects  of  double  passage  of 
waves  in  randomly  inhomogeneous  media,”  Sov  Phys  Usp 
25. 494-508(1982). 

4.  E,  Jakeman,  “Enhanced  backscattenng  through  a  deep  ran¬ 
dom  phase  screen,"  J.  Opt.  Soc.  Am.  A  5,  1638-1648 1 1988; 

5.  P.  R.  Tapster.  A.  R.  Weeks  and  E  Jakeman.  "Observation  of 
backscattenng  enhancement  through  atmospheric  phase 
screens,”  J.  Opt.  Soc.  Am  A  6,  517-522 1 1989). 

6.  R.  Mazar  and  A.  Bronshtein,  "Double  passage  analysis  in 
random  media  using  two-scale  random  propagators.”  Waves 
Random  Media  1, 341-362 1 1991 1. 

7.  E.  Jakeman.  P.  R.  Tapster,  and  A,  R.  Weeks,  "Enhanced 
backscattering  through  a  deep  random  phase  screen,"  J  Phys, 
0  21,32-36,1988). 

8.  G.  Welch  and  R.  Phillips.  "Simulation  of  enhanced  backscatter¬ 
ing  by  a  phase  screen,”  J.  Opt.  Soc.  Am.  A  7,  578-584  1 1990). 

9  B.  S.  Agrovskii,  A.  N.  Bogaturov.  A.  S.  Gurvich,  S.  V.  Kireev, 
and  V.  A.  Myakinin,  "Enhanced  backscattering  from  a  plane 
mirror  viewed  through  a  turbulent  phase  screen.”  J.  Opt.  Soc. 
Am.  A  8, 1142-1147(1991). 

10.  G,  J,  BaJmer,  D  L.  Jordan,  P,  R.  Tapster,  M.  J.  Kent,  and  E 
Jakeman,  “Double  passage  wave  propagation  effects  at  infra¬ 
red  and  visible  frequencies,”  Opt.  Commun.  88,  6-12 1 1992). 

11.  J.  W.  Goodman,  "Statistical  properties  of  speckle  patterns,”  in 
Laser  Speckle  and  Related  Phenomena,  2nd.  ed..  J.  C.  Dainty, 
ed.  (Springer-Verlag,  Berlin,  1984),  pp.  9-75. 

12.  D.  L6ger,  E.  Mathieu,  and  J.  C,  Perrin,  “Optical  surface 
roughness  determination  using  speckle  correlation  tech¬ 
nique,”  Appl.  Opt.  14, 872-877  ( 1975). 

13.  D.  L6ger  and  J.  C.  Perrin,  "Keai-time  measurement  of  surface 
roughness  by  correlation  of  speckle  patterns,”  J.  Opt.  Soc.  Am. 
66, 1210-1217(1976). 

14  S.  Feng,  C.  Kane,  P.  A.  Lee,  and  A.  D.  Stone,  “Correlations  and 
fluctuations  of  coherent  wave  transmission  through  disor¬ 
dered  media,”  Phys.  Rev.  Lett.  61, 834-837  (1988). 

15.  R.  Berkovits,  M.  Kaveh,  and  S.  Feng,  “Memory  effect  of  waves 
in  disordered  systems;  a  real-space  approach,”  Phys.  Rev.  B 
40, 737-740(1989). 

16.  S.  Feng,  “Novel  correlations  and  fluctuations  in  speckle  pat¬ 
terns,”  in  Scattering  and  Localization  of  Classical  Wooes  in 
Random  Media,  P.  Sheng,  ed.  (World  Scientific,  Singapore, 
1990),  pp.  179-206. 

17.  I.  Freund  and  M,  Rosenbluh,  “Time  reversal  symmetry  of 


2742  APPLIED  OPTICS  /  Vol.  32,  No.  15  /  20  May  1993 


multiply  scattered  speckle  patterns,”  Opt.  Commun  82,  362- 
3C-  ,19911. 

18.  C.  J.  R.  Sheppard  and  T.  Wilson.  "Image  formation  in  scan¬ 
ning  microscopes  with  partially  coherent  source  and  detector,  ' 
Opt.  Acta  25, 315-325  ( 1978). 

19.  E.  Jakeman  and  J.  G.  McWhirter,  "Correlation  function 
dependence  n.  -.ae  scintillation  behind  a  deep  random  phase 
screen."  ■'  ‘’hys.  A  10,  1599-1643 1 19771. 

20.  I  S  i,e>a,  "On  a  moment  theorem  for  complex  Gaussian 
p'-'j„  ^oes,”  IRE  Trans.  Inf  Theory  IT-8,  194-195  1 1962). 

21  M.  P.  Van  Albada  and  A.  Lagendijk,  "Observation  of  weak 
localization  of  light  in  a  random  medium,"  Phys.  Rev.  Lett  55, 
2692-269511985). 

22.  P.  E.  Wolf  and  G.  Maret,  "Weak  localization  and  coherent 


backscattering  of  photons  in  disordered  media,”  Phys  Rev 
Lett.  55,  2696-2699  1985 : 

23  A.  Ishimaru,  "Experimental  and  theoretical  studies  of  en¬ 
hanced  backscattering  from  scatterers  and  rough  surfaces,  in 
Scattering  in  Volumes  and  Surfaces,  M  .Nieto- Vesperinas  and 
J  C,  Dainty,  eds  iNorth-Holiand.  Amsterdam,  1990 1.  pp 
1-15 

24  1  Freund.  M  Rosenbluh,  and  R  Berkovits.  'Geometric  scaling 
of  the  optical  memory  effect  in  coherent-wave  propagation 
through  random  media,"  Phys  Rev  B  39,  12403-12406 
<l989i 

25.  P  Beckmann  and  A.  Spizzichino.  The  Scattering  of  Electroniag- 
netic  Waves  from  Rough  Surfaces  Pergamon,  Oxford,  1963  . 
pp  85-87 


1  Accesion  Fo7 

V- 

NTIS 

CRA&I 

OTIC 

TAB 

□ 

Unannounced 

□ 

Justification 

Bv 

Distribution  | 

Availability  Codes 

Avail  and  /  or  1 

Oist 

Spe 

cial 

A-i 

dsO 

jynC  QUALITY  inspected  2 


20  May  1993  Vol.  32,  No.  15  APPLIED  OPTICS  2743 


